AD-A185  718 


UTIC  fill  CUPJI 


ncorjmTti  nopggt  (Mt-er-i 


Asniraoiic  toKtixa  of  fqut-t  isiBmEs 

MEm  EBSESDH  ARUOKCKHS 


J^ 


G0GD^S  1l%  WCKl^ 

Divisicn  of  Biostatistics 
MeBocial  Sloan-iBsttering  Cianaer  Genter 
Nev  Yodc,  New  Yodc  10021 
Ocbcber,  1987 


OTIC 

ELECTE 
OCT  2  3 1987L 


C^i 


Prepeared  under  ocntract  NO.  N00014-85-K-0485,  NR  150-536 
with  the  Personnel  eoid  TTedning  Research  Progreons 
Psychological  Sciences  Division 
Office  of  Naval  Research 


I  A(]proved  for  public  release;  distributicai  unlimited. 

ReproductiCTi  in  whole  or  in  peart  is  permitted  for 
.  any  purpose  of  the  Uhited  States  Government. 


8  7  15 


unclassified 


m 


REPORT  DOCUMENTATION  PAGE 

la.  REPORT  SECURITY  CLASSIFICATION 

unclassified 

1b  RESTRICTIVE  MARKINGS 

2a.  SECURITY  CLASSIFICATION  AUTHORITY 

unclassified 

3.  DISTRIBUTION /AVAILABILITY  OF  REPORT 

Approval  for  public  release;  distribution 
unlimited 

2b.  DECLASSIFICATION /DOWNGRADING  SCHEDULE 

4.  PERFORMING  ORGANIZATION  REPORT  NUM8ER(S) 

S.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

6a.  NAME  OF  PERFORMING  ORGANIZATION 

George  Y.  Wong 

6b.  OFFICE  SYMBOL 
(If  applicable) 

7a  NAME  OF  MONITORING  ORGANIZATION 

Office  of  Naval  Research 

6c^ORESS  (pty.  State,  and  ZIP  Code) 

Memorial  Sloan-Kettering  Cancer  Center 

1275  York  Avenue 

New  York,  New  York  10021 

7b.  ADDRESS  (C/ty,  State,  and  ZIP  Code) 

800  North  Quincy  Street 

Arlington,  VA  22217 

8a.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 

Personnel  &  Training  Research 

8b.  OFFICE  SYMBOL 
(If  applicable) 

9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

NOOO  14-85-K-0485 

8c ^DRESS  (Oty.  State,  and  ZIP  Codt) 

Office  of  Naval  Research,  Code  442 
800  North  Quincy  Street 
Arlington,  VA  22217 


11.  TITLE  Onclude  Security  Classification) 


10.  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM  PROJECT 

ELEMENT  NO.  NO. 

NR150>536 


Asymptotic  Normality  of  Poly-t  Densities  With  Bayesian  Applications 


12.  PERSONAL  AUTHOR{S)  ..  , 

George  Y.  Wong 


5.  PAGE  (COUNT 


16.  SUPPLEMENTARY  NOTATION 


COSATI  CODES 


riuctnn  ) 


18  SUBJECT  TERMS  {Continue  on  reverse  if  rtecessary  and  kfentify  by  block  number) 

I  CC^VJ.ln  /  - -  y- 


^  *  A  \ _  • _ ■  J 


19.  ABSTRACT  {Continue  on  reverse  if  necessary  and  identify  by  l^ock  number)  ^  _  ^  _  _ 

^  A  poly-^  density  Is  a  density  ^itajdi  j^p^rad^rtioned  to  a  product  of  ^  least  two 
t-lite  factors,  eadti  of  vAiich  is  of  -ttke  wiere  d  is  a 

positive  nuDber,  is  an  arbitrary  loraticn  vector  and  D  is  a  syninefbric  semi-positive 
definite  scale  matrix.  In  general,  ^  is  a  function  of  d.  Such  a  dsisity  arises,  for 
exaoople,  in  the  Bayesian  analysis  of  a  linear  model  with  a  normal  error  term,  independ^ 
normal  priors  on  the  linear  parameters  and  inverted-ganma  priors  on  the  variance 
ccnponents.  A  theorem  about  the  asynptotic  normality  of  the  density  as  a  sidoset  of  the 
individuEd  d's  tend  to  infinity  is  proved  under  very  general  conditions.  A  corollary 
^lecificedly  related  to  the  Bayesian  linear  model  is  also  given.  Detailed  results  are 
illustrated  in  the  familiar  Bayesian  multiple  linear  regressicn  model  with  two  variance 
conponarts.  The  Tiao-Zellner  e^qsansion  for  approximating  the  particular  poly-t  form 
involving  two  proper  multivariate  t  factors  is  extended  to  the  case  of  two  arbitrary 


I O I  •  "ST! 


20.  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT 

fiOuNCLASSIFIED/UNLIMITED  □  SAME  AS  RPT.  □  DTIC  USERS 

21.  ABSTRACT  SECURITY  CLASSIFICATION 

unclassified 

22a 

NAME  OF  RESPONSIBLE  INDIVIDUAL 

Dr.  Charles  Davis 

|2b  TELEPHONE  (Include  Area  Code) 

_ ^02/696-4046 _ 

22e.  OFFICE  SYMBOL 

DO  FORM  1473, 84  MAR 


83  APR  edition  may  be  used  until  exhausted. 
All  other  editions  are  obsolete. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


DTIC  TAB 

Unannounced 

Justification 


□ 

□ 


ASYMPTOTIC  NORMALITY  OF  POLY-T  DENSITIES 
WITH  BAYESIAN  APPLICATIONS 


GEORGE  Y.  WONG 


DIVISION  OF  BIOSTATISTICS 
MEMORIAL  SLOAN-KETTERING  CANCER  CENTER 
NEW  YORK.  NY  10021 


By . . . . 

Dii.t:‘ibution/ 


Av.^iiabitity  Codes 

'  Avo'!  and/or 
Gptyrial 


Key  Words  and  Phrases:  poly-t  density;  asymptotic  normality;  generalized  Tiao- 

Zellner  expansion;  Bayesian  linear  model;  posterior 
inference  of  linear  parameters. 


ABSTRACT 

A  poly-t  density  is  a  density  which  is  proportional  to  a  product  of  at 
least  two  t-like  factors,  each  of  which  is  of  the  form  {1  ♦  (x  -  ij)^M 
where  d  is  a  positive  number,  v  is  an  arbitrary  location  vector  and  M  is  a 
symmetric  semi-positive  definite  scale  matrix.  In  general,  M  is  a  function  of 
d.  Such  a  density  arises,  for  example,  in  the  Bayesian  analysis  of  a  linear 
model  with  a  normal  error  term,  independent  normal  priors  on  the  linear 
parameters  and  inverted-gamma  priors  on  the  variance  components.  A  theorem 
about  the  asymptotic  normality  of  the  density  as  a  subset  of  the  individual  d's 
tend  to  infinity  is  proved  under  very  general  conditions.  A  corollary 
specifically  related  to  the  Bayesian  linear  model  is  also  given.  Detailed 
results  are  illustrated  in  the  familiar  Bayesian  multiple  linear  regression 
model  with  two  variance  components.  The  Tiao-Zellner  expansion  for 
approximating  the  particular  poly-t  form  involving  two  proper  multivariate  t 
factors  is  extended  to  the  case  of  two  arbitrary  t-like  factors. 


1.  INTRODUCTION 


A  p-dimenslonal  random  vector  X  is  said  to  have  a  po1y-t  distribution  if 
its  density  is  proportional  to  a  product  of  L>2  "t-like"  factors,  or 

f(x|dj . djo  n  {l>(x-ji^)\(d,^)  (x-]i^)}.’‘^k'2  jjj 

I^I 

where  d|^>0,  is  a  pxl  location  vector  and  M|^(d|j)  is  a  pxp  symmetric  semi¬ 
positive  definite  scale  matrix,  k=l,...,L.  To  ensure  that  the  right-hand  side 
of  (1),  which  we  denote  by  g(xldj,...,dj^),  is  normalizable,  we  must  require 
that  dj+...+d^>p  and  Hj(dj)+. ..♦M|^(dj^)  be  positive  definite  (Dickey  1986). 

The  normalizing  constant  fg(x.|dj,...,d|^)dx,  however,  cannot  be  expressed  in  a 
simple  closed  form.  When  L  is  smaller  than  p,  Dickey  showed  that  the 
normalizing  constant  can  be  expressed  in  an  (L-l)-fold  integral. 

A  useful  form  of  density  (1)  is  obtained  by  letting 

V'’k  *  "k-  fikt^k)  *  Sk''’k  (2) 

for  some  V|^>0,  mj^  such  that  mj^>0,  and  Hj^  symmetric  semi-positive  definite. 
When  m|^=  p  and  is  positive  definite,  then  the  kth  t-like  factor  is 
proportional  to  a  proper  multivariate  t  density  with  degrees  of  freedom. 

Such  a  poly-t  density  form  plays  an  important  role  in  the  Bayesian  analysis  of 
a  general  linear  model  with  a  normal  error  term.  The  Bayesian  approach 
considered  here  is  described  as  follows:  (1)  A  vague  prior  is  imposed  on  a 
subset  of  the  linear  parameters;  (2)  The  rest  of  the  parameters  are 
partitioned  into  different  subsets,  and  an  exchangeable  normal  prior  with  mean 
0  and  an  unknown  variance  component  is  imposed  on  each  subset  independently;  and 
(3)  The  prior  variances  and  the  error  variance  are  given  independent  inverted 
gamma  distributions.  It  can  be  shown  that  the  posterior  density  of  the  linear 
parameters  is  of  the  above  poly-t  density  form.  Details  about  poly-t  posterior 


distributions  can  be  found  In  Dickey  (1974),  Oreze  (1977),  Rajagopalan  and 
Broemeling  (1983),  Broemellng  and  Abdullah  (1984),  and  Broemellng  (1985). 

In  this  article,  we  are  concerned  with  the  asymptotic  normality  of  a  poly- 
t  density  as  some  or  all  of  dj,...,dj^  become  large.  The  asymptotic  normal 
density  can  be  used  to  approximate  the  poly*t  density  directly.  Also,  It  can 
serve  as  an  Important  sampling  distribution  for  the  Monte  Carlo  evaluation  of 
any  probability  statement  about  X,  or  any  moment  of  X.  Despite  Its  useful Iness 
In  Bayesian  linear  modeling,  the  asymptotic  normality  of  a  poly-t  density  has 
never  been  formally  established.  In  Section  2,  we  prove  the  normality  result 
for  the  po1y-t  density  (1)  under  fairly  general  conditions.  As  will  be  seen 
there,  the  proof  Is  much  harder  than  that  In  the  case  of  a  single  multivariate 
t  density;  the  difficulty  is  mainly  due  to  the  fact  that  the  normalizing 
constant  of  a  poly-t  density  Is  In  terms  of  an  Intergral  expression. 

In  Section  3,  we  illustrate  the  asymptotic  results  on  a  Bayesian  linear 
model  with  two  variance  components.  There  we  show  that  the  appropriate 
asymptotic  normal  density  may  be  made  the  leading  term  In  a  generalization  of 
the  Tiao-Zellner  expansion  (1964)  to  approximate  the  posterior  density  of  the 
linear  parameters. 

2. MAIN  RESULTS 

We  show  that  under  some  fairly  general  conditions,  the  poly-t  density  (1) 
converges  to  a  proper  multivariate  noraml  density.  It  follows  from  a  result  of 
Scheffe  (1947)  that  the  poly-t  variable  actually  converges  in  distribution  to 
the  multivariate  normal  variable  with  the  limiting  normal  density.  We  also 
give  a  corollary  on  a  special  form  of  (1)  which  is  particularly  useful  in 
Bayesian  linear  modeling.  For  convenience,  we  write  A>0  and  A>0  to  mean  that 
the  symmetric  matrix  A  is  semi-positive  definite  and  positive  definite, 
respectively.  Also,  we  denote  the  density  of  a  multivariate  normal  density 


with  the  mean  p  and  covariance  matrix  £  by  ). 

Theorem  For  the  general  poly-t  density  (1),  assume  that 

(1)  converges  to  M|^>0  as  •»  k*!,...,  L,  and  M  *  Mj 
and 

(2)  there  exist  semi-positive  definite  matrices  Bj,...,  such  that 

B|^>0  and  ■“A  (dj^)  -  ^>0  for  all  dj^  sufficiently  large. 

Then  as  dj,...,dj^  all  tend  to  infinity,  f(x|dj,. . .d^)  will  tend  to  the  limiting 
density  ff|(x|w,H"^) ,  where 

Proof.  Without  any  loss  of  generality,  we  prove  the  theorem  for  the  case 
L*2.  To  show  that  f{xidj,d2)  converges  to  as  dj,  d^-^",  it 

suffices  to  show  that  (i)  g(x|dj,d2)  tends  to  a  limit  g(x)  proportional  to 
F|^(xlil.M  ^)»  and  (ii)  fg(x|dj,d2) -►  fg(x)dx,  as  dj,  d2-^». 

To  prove  (i),  we  make  use  of  assumption  (1)  and  obtain 

2  H  n 

g(x!di.d2)  ♦  Q,^/d^  ♦  o(l/d^)}'VS 

where  Qj^  *  (x-ii,^)^M|^  (-"Hfe)*  k  =  1,  2.  It  follows  that  g(xldj,d2)  tends  to 
g(x)  =  exp  {  -  x  exp  {  -  (x-v)^M(x-ij)/2} , 

which  is  proportional  to  f|^(x!]iiM”^). 

To  prove  (ii),  we  show  that  for  dj  and  d2  sufficiently  large,  g(xjdj,d2) 
is  dominated  by  an  integrable  function  h(xjdj,d2),  which  monotonically 
decreases  to  an  integrable  function  h{x);  moreover,  fh{x [dj,d2)d£  also 
decreases  monotonically  to  fh(x)dx.  Then  using  a  generalization  of  the 
Lebesque  Dominated  Convergence  Theorem  (see,  for  example,  Royden,  p.89),  we 
conclude  that  the  limit  of  fg(x|dj,d2)dx  is  fg(x)dx. 

We  now  make  use  of  assumption  (2)  to  prove  the  above  claims.  First,  it  is 
an  immediate  consequence  of  the  assumption  that  for  sufficiently  large  dj  and 

g(x|dj,d2)  i  h(x|dj.d2)  *  n  {1  ♦  (x-y,^)\(x-Hk)/dk}'‘^k^2^ 


Also,  it  can  be  directly  verified  that  as  dj,d2-»-,  h(x|dj,d2)  monotonically 
decreases  to  an  integrable  function 

h(x)  -  exp  {  - 

X  exp  {  -  (X-  v)^(Bj+B2)  (x-v)/2}, 

where  v  ■  (Bj*B2)”^(Bjjjj+B2ii2) .  To  establish  the  integrabil ity  of  h(j(ldjd2), 
let  m  »  min(dpd2)>  Then  the  raonotonicity  of  the  function  implies  that 

2 

h(xdj,d2)  <^^{1  ♦ 

<{1  ♦  (x-v)^(Bj*B2)(x-v)/m)‘'"'^ 

For  sufficiently  large  dj  and  d2,  the  t-like  factor  on  the  right-hand  side, 
denoted  by  h^(x),  is  proportional  to  a  multivariate  t  density  with  m-p  degrees 
of  freedom,  mean  v  and  scale  matrix  ^(x|<ii*^2^  integrable. 

Let  m  be  fixed.  Define  ♦(x|dj,d2)  =  -  h(x|dj,d2),  for  d^,  d2>m. 

Applying  the  Monotone  Convergence  Theorem  to  this  sequence  of  functions 
together  with  the  fact  that  hj|j(x}  is  integrable,  we  conclude  that 
fh(x|dj,d2)dx  tends  to  fh(x)dx).  This  completes  the  proof. 

We  point  out  that  if  only  the  elements  of  a  subset  of  dj,...,dj^  tend  to 
infinity,  then  the  theorem  will  apply  to  the  product  of  the  corresponding  t- 
like  factors  while  the  remaining  factors  will  stay  put. 

We  now  consider  a  special  form  of  the  poly-t  density  (1)  with  scale  matrix 

W  ‘  "k'Okifik  (3) 

for  some  positive  function  Wj^  and  semi-positive  definite  matrix  The 
following  corollary  concerns  the  asymptotic  behavior  of  such  a  poly-t  density. 

Corollary  For  the  particular  poly-t  density  with  M|^(d|^)  of  the  form  (3), 
k»l,...,L,  assume  that  C|^>0  as  dj^-*-  •,  k=l,...,L  and  M  =  Mj+. .  .♦Mj^>0, 


4 


Mhere  Mj^  ■  'A  the  limit  of  (d,^).  Then  f(iidj,...,dj^)  as 

dj,. . .  ,dj^-^». 

Proof  We  need  to  check  that  the  two  asusmptions  of  the  theorem  hold  for 
this  particular  form  of  poly-t  density.  Assumption  (1)  is  obviously  satisfied. 
For  assumption  (2),  we  take  ^  appropiately  chosen  positive 

number  Clearly,  We  must  choose  such  that 

dkM|^(dk)-Bj^  =  {d|jW|^(dj^)  -  *^|5Cjj}N^>0  for  dj^  sufficiently  large.  Since  d|^W|^{d|^)-^ 
C|^  as  d|^-^«,  it  follows  that  for  dj^  large  enough,  d|^W|^(d|^)  ^  Cj^  -  n|^>0  where 
is  an  arbitrary  positive  number.  Therefore,  we  can  choose  Ej^  =  (C|j-n|^)/C|^. 
This  completes  the  proof. 

An  example  of  a  poly-t  density  with  scale  matrices  of  the  specific  form 
(3)  is  when  dj^  and  M|j(d|^)  are  given  by  (2).  Here  we  may  choose  W|^(d|^)  = 
mj^),  and  Nj^  “  ^  matrix  in  (2)  is  the  limit  of  d|^M^^(d|^) 

since  Cj^  =1. 

3.  APPLICATION  TO  BAYESIAN  LINEAR  MODELING 

The  results  of  the  previous  section  are  now  applied  to  a  multiple 
regression  model  with  an  exchangeable  normal  prior  on  the  regression 
coefficients  and  independent  inverted  gamma  priors  on  the  variance  components. 
The  sampling  model  is  represented  by 

1  1  *£  » 

e  -v.  N(0,o%), 

where  ^  is  an  nxl  vector  of  observations,  X  is  an  nxp  fixed  design  matrix  in 
correlation  form,  and  I^  is  the  nxn  identity  matrix.  The  prior  specification 
is  given  by 


i  N(ylp.  o| 

vX/o^  'V.  xj.  xj  , 

0 

Where  Ip  is  a  px!  vector  of  ones,  and  the  x’s  and  v's  are  hyperparameters  that 
determine  the  prior  means  and  variances  of  the  variance  components.  For 
example  E(o^)  =  vX/(v-2),  and  V(o^)  *  2vV/(v-2)^  (v-4). 

Therefore,  a  large  value  of  reflects  a  strong  prior  certainty  the  X  is  a 
correct  guess  of  o  .  Additionally,  we  assume  that  prior  knowledge  about  y  is 
vague. 

In  the  special  case  of  the  one-way  random  effects  model  with  a  block 
diagonal  matrix  X  =  <ii3g(^,. . . ,1  ),  where  n.  is  the  sample  size  of  the  ith 

IP 

group.  Hill  (1965,  1977,  1980)  presented  exact  and  approximate  posterior 
2  2 

inference  of  and  o  .  His  results  can  be  used  to  obtain  the  posterior 
moments  of  the  group  means  8.  Lindley  and  Smith  (1972)  estimated  the 
parameters  of  the  above  general  model  using  a  computationally  simple  posterior 
joint  modal  approach.  None  of  these  authors,  however,  discussed  the 
approximation  of  the  posterior  density  of  8,  which  we  denote  by  f(8l^).  In 
the  following,  we  show  how  the  results  of  Section  2  may  be  used  to  approximate 
this  density. 

From  Lindley  and  smith,  the  posterior  density  of  8  is 

^(ilK)  o  (l*iWvgXg)"^''8  X  {l*(l-i)Vx(8-8)/(vX  +ns^)}’^'^*"^'^, 

where  H  =  lp"(l/p)lplp^,  8  is  a  least  squares  estimate  of  8  and  s^  is  the 
residual  mean  square.  The  posterior  density  of  8  is  a  poly-t  density  of  the 
special  form  (2)  with 

2 

vi  *  Vg»  V2  *  v  ns  /x, 
n>j  *  p-1,  m2  =  n  -  ns^/X, 

Mj  »  H/Xg,  *  xV^- 


2  7 

When  there  is  prior  near  certainty  concerning  and  o  so  that  both  v. 

P  P 

and  V  are  large,  a  simple  application  of  the  corollary  shows  that  f(£i^)  is 
approximately  normal  with  mean  and  covariance  matrix  z  given  by 


u  »  (H/Xg  ♦ 

Z  =  (H/Xg  ♦  X^XA)'^ 

When  there  is  a  considerable  amount  of  uncertainty  concerning  the  prior 
value  X  (orx)  so  that  v  (orv)  is  not  assigned  a  large  value,  the  normal 

M  P 

approximation  may  be  improved  upon  by  using  a  generalization  of  the  Tiao- 
Zellner  expansion.  An  alternative  was  proposed  by  Dickey  (1967)  using  an 
expansion  based  on  Appell's  polynomials  (1880)  in  the  univariate  case  (p=l); 
however,  the  multivariate  version  has  not  been  studied.  The  original  Tiao- 
Zellner  expansion  is  designed  for  a  special  case  of  a  poly-t  density  with  L=2 
proper  multivariate  t  factors  with  mean  ji|^,  scale  matrix  and  V|^  degrees  of 
freedom,  k*l,2.  In  the  expansion,  the  poly-t  density  is  approximated  by  a 
density  in  the  form  of  fm(x{ii,M"^) ,  where  ji  =  M"^(Mjiij+M2ii2)  M  =  Hj  ♦ 

^2,  multiplied  by  a  double  infinite  series  in  the  inverse  powers  of  vj  andv^* 
When  most  of  the  probability  mass  is  concentrated  in  the  intersection  of  two 
ellipsoids  {x|  (x-iJ|^)^M|^(x-jij^)  i  vj^},  k*l,2,  the  expansion  will  yield 
reasonable  approximation  with  only  a  few  terms  in  the  series  expansion. 

The  Tiao-Zellner  expansion  may  be  generalized  to  the  case  of  a  poly-t 
density  of  the  form  (2).  Lindley  (1971)  mentioned  the  use  of  this  expansion  in 
the  case  of  the  one-way  model.  However,  he  did  not  specify  the  limiting  normal 
density.  Also,  the  degrees  of  freedom  vj  and  V2  were  not  correctly  stated.  We 
now  sketch  the  essential  steps  in  the  generalization  using  f(B|^)  in  (4)  as  an 
example. 

Let  V|^  and  mj^  be  as  defined  in  (5),  Qj  =  and  Q2  =  (6- 

T  T 

i)  X  X(8-6)/X.  Following  Tiao  and  Zellner,  we  obtain  for  the  first  factor 


.  exp(-Qj/2)  exp[Qj/2-{(Vj>mj)/2}  log 

This  is  of  the  same  form  as  that  obtained  In  (3.2)  of  Tiao  and  Zellner  except 
that  we  replace  p  there  by  m^.  Therefore,  the  functional  forms  of  the 
polynomials  p^  *  P^CQl)  in  the  expansion 

00  4 

=  exp(-Q,/2)  I  P4  vj' 
j=0 

are  of  the  same  forms  as  those  given  in  Tiao  and  Zellner  except  that  the 
constant  p  is  replaced  by  mj,  (for  instance,  pj  =(Qj-mjQj)/4) .  Similarly,  for 
the  second  factor. 


(l+Q./v  )'^V'"2^^^  =exp(-Q./2)  z  q 

j=0  j  2 

where  the  constant  p  which  appears  in  the  polynomials  q.  =  q.(Q2)  in  Tiao  and 

J  J  t 

Zellner  is  now  replaced  by  m2. 

Following  the  rest  of  the  derivation  of  Tiao  and  zellner,  we  obtain 
the  generalized  expansion 


i-1 


CD  es 


“1.  -0 


1=0  1=0  ^ 


j=0 

where  the  leading  normal  term  is  the  one  guaranteed  by  the  corollary,  and  the 
d.j. 's  are  polynomials  in  p^  and  q^  given  in  (3.9)  of  Tiao  and  Zellner.  The 
coefficients  which  appear  in  the  polynomial  d..,  namely  b  =  E{p  (Q,)q  (Q,)} 
for  r<i,  s£j,  also  have  to  be  modified  accordingly.  The  modification 
essentially  consists  of  the  following  steps:  (1)  Change  the  constant  p  in 
p^(Qj)  given  in  Tiao  and  Zellner  to  mj;  (2)  Change  the  constant  p  in  Q^(Q2) 
given  in  Tiao  and  Zellner  to  m2,  and  (3)  Use  the  bivariate  moment-cumulant 
inversion  formulas  as  given  by  Cook  (1951)  to  express  b^^  in  terms  of  the  mixed 
cumulants  of  Qj  and  Q2.  Expressions  for  these  cumulants  are  derived  in  Tiao 
and  Zellner.  Following  these  four  steps,  the  coefficients  b^^  can  be  evaluated 


in  a  straightforward  manner.  We  omit  the  details  here. 

In  the  case  of  a  product  of  two  t-like  factors  of  the  form  (2),  it  is 
also  feasible  to  approximate  the  density  by  direct  one-dimensional  numerical 
integration.  Dickey  (1968)  showed  that  integration  can  be  facilitated  if  the 
scale  matrices  are  simultaneously  diagonalized  first  (more  details  are 
presented  in  Box  and  Tiao  (1973,  chapter  9)).  The  generalized  Tiao-Zellner 
expansion  offers  a  simple  approximation  alternative  in  terms  of  familiar 
multivariate  normal  calculations.  When  there  are  many  such  t-like  factors, 
such  as  in  the  posterior  analysis  of  a  general  mixed  linear  model  with  many 
variance  components,  it  will  be  virtually  impossible  to  approximate  the 
density  by  high-dimensional  numerical  integration.  The  general ized 
expansion,  however,  can  be  extended  to  the  case  of  many  t-like  factors  in  a 
straightforward  manner.  Obviously,  the  computational  complexity  will 
increase  rapidly  as  the  number  of  t-like  factors  increases. 
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